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The existence of a finite electric dipole moment (EDM) d a of the closed-shell atom Yb implies 
parity and time reversal violations involving the nuclear sector. An important effect which can 

S contribute to the Yb EDM is the tensor-pseudotensor electron-nucleus interaction characterized by 
the coupling constant Ct- Within the Standard Model (SM) of particle physics Ct = 0, as this form 
I *; ■ of interaction is not allowed. If a finite d a of Yb is observed in experiments, then an estimate of 

Ct can be obtained by combining with the theoretical calculations. A non-zero Ct implies physics 
^3 , beyond the standard model. In this paper we present the result of our ab initio calculation of the 

O ' EDM Yb using different many-body methods. 
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^ ■ I. INTRODUCTION 

Discrete symmetry violations in atoms are important phenomena to probe for physics beyond the Standard 
. Model(SM) of the particle physics. The electric dipole moment(EDM) of an atom, which is a signature of the 
simultaneous parity and time-reversal (P-T) symmetry violations is one such. There are several possible sources of 
P-T violation effects within an atom, the tensor pseudo-tensor electron nucleus interaction is an example which is 
semi-leptonic in nature. The atomic Yb, which is a closed-shell atom having a Z = 70 is a very good candidate for 
, atomic EDM experiments to probe for the nuclear sector effects. An additional advantage of using a rare Earth atom 
like Yb is the closely spaced energy levels. The experimental results when compared with the theoretical results can 
yield signatures of physics beyond the SM. 

Among the closed-shell atoms, the atom EDM of Xe jl| and Hg || have been measured. Though the Hg experiment 
has set the record of being the most sensitive spectroscopy ever done, the result obtained is null and sets an upper 
bound to atomic Hg EDM as 2.1 x 10~ 28 ecm. Still, it has set bounds on the parameters in particle physics models [|J. 
Improving the accuracy further is an important challenge in atomic EDM experiments. Atomic Yb offer a possibility 
^*"> of achieving this by using the techniques of laser cooling and trapping. It is also desirable that the EDM measurements 
• be done in other closed-shell atoms to verify the physical effects observed. 

Atomic Yb also has the advantage of being relatively simple in level structure among the rare Earth elements. It has 
^ • been studied theoretically using a variety of atomic many-body methods and there are ongoing experimental studies. 
' The exc itation energies has been studied using relativistic coupled-cluster method ; the hyperfine structure constants 
and electric-dipole transition properties have been studied using multireference relativistic many-body perturbation 
theory ^|,[5| ; and the oscillator strength of some of the important transitions were investigated using multi-configuration 
relativistic Hartree-Fock ||. The lifetime of some of the crucial levels have been measured experimentally and 
atomic Yb has been laser cooled and trapped |BS|. It is the candidate for the atomic EDM experiments using the 
methods of laser cooling and trapping of atoms which are in progress Jl0|,[ll} • 

The P-T violating effects which can be probed using atomic Yb are the tensor pseudotensor (T-PT) electron-nucleus 
interaction |l^] and the Schiff moment p3[ . An observation of a finite EDM of atomic Yb could mean nonzero Ct-i 
which is a signature of physics beyond the SM as Ct is zero within SM. Hence the theories which allow T-PT electron- 
nucleus interactions would play a role in understanding nature. The parameter Ct is extracted from the experimental 
result of atomic EDM by combining with the theoretical calculations. Accurate atomic theory calculations are needed 
to obtain precise Ct .The atomic EDM can also put bounds on the parameters of alternative models in particle physics. 
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The Schiff moment which arises due to either nucleon EDM or P-T violating interactions between the nucleons can also 
contribute to the Yb EDM and can be used to extract the nuclear P-T violation parameters. We have calculated the 
contribution to atomic Yb EDM from both these effects but in this paper we will present the result of our calculations 
for T-PT contribution alone, in a later paper we will present the result of our Schiff moment Computation. 

The calculation of Yb EDM within a limited set of configuration state functions (CSFs) using diagonalization and 
Bloch equation based many-body perturbation theory ( MBPT) methods was reported in an earlier paper JlJ|. For 
convenient of reference this paper is called as paper-I hereafter. As a followup, in this paper we present results of 
calculations using much larger CSF space. In addition, we discuss the drawbacks common to the diagonalization 
and Bloch equation based perturbation methods, improvements are discussed and results of using these methods are 
presented. In this paper we present the method of the calculation. 



II. BLOCH EQUATION BASED MANY-BODY PERTURBATION THEORY 



This section is a brief overview of the interaction used in the calculation and Bloch equation based MBPT. A 
detailed description of which are given in the paper-I. The effective T-PT electron-nucleus interaction Hamiltonian 
obtained by treating the nucleus nonrelativistically is 



(2.1) 



where Ct is the T-PT electron- nucleus coupling constant, Gf is the Fermi coupling constant, / is the nuclear spin, /3 
and a are the Dirac matrices and pn(t) is the nuclear density. This interaction Hamiltonian is effective only within 
the nuclear region, where pAr(r) is non-zero. And the dependence on / implies that it is observable only in odd 
isotopes of Yb which has nonzero /. We use Fermi nuclear density model for our calculations. 

The Bloch equation of an atom with unperturbed Hamiltonian Hq and the residual coulomb interaction V es as the 
perturbation is 



[f2 cs , Hq] — Vcs^cs — XcsPV cs Q c 



(2.2) 



Where Q os is the wave-operator, Xcs = ^es — 1 is the correlation operator, P and Q are the projection operators 
of the model and complementary space. The first and second terms on the right hand side are the principal and 
renormalization terms respectively. Using Epstein-Nesbet(E-N) partitioning Hq and V es can be defined in terms of a 
set of CSFs as 



#0 


= E 

i 




-^atom 


*> 


*■> 


<*, 


+£< 

i 




^atom 


*,) 




(*< 


v cs 


=£ 


;*< 


^atom 




*,) 




+E 




■"atom 









where H atom is the Dirac-Coulomb atomic Hamiltonian, and and {l^i)} are the even and odd parity CSF 

spaces respectively. The Dirac-Coulomb Hamiltonian of an atom with N electrons in atomic units ( e = 1, h = 1, 
and m e = 1 ) is 
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(2.3) 



where a,; is the Dirac matrix, pi is the momentum of the electron and V^uc is the nuclear potential. Introduce the 
P-T violating interaction Hamiltonian H FTV as perturbation and define the total perturbation Hamiltonian V as 



V = V es + H p 



(2.4) 



The interaction Hamiltonian H PTV can also be expressed in terms of CSFs as in Hq and V cs . Treating V as the 

s the total wave operator 

(2.5) 



perturbation and defining fi(edm) as the total wave operator 

fi(edm) = tt os + f2 GS ,edm, 
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where Q es remains the same and fi C s.cdm is the wave-operator which has one order of H rTV and all possible orders of 
14s before and after H PTV . The H rTV is treated to first order since it scales as Gf, which is very small. The expression 
of the atomic EDM d a is 
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(2.6) 



where D c s = ilJ s Z3r2 cs cdm + iY cdm -Df2 cs is the effective atomic EDM operator. It is the dipole operator dressed 
with the all order residual coulomb interaction and one order of H PTV arranged in all possible sequence. The Bloch 
equation based MBPT has been used for the calculation with a large CSFs space as it is more efficient in terms of 
execution time unlike the direct matrix diagonalization approaches. 



III. SIZE CONSISTENT THEORY IN CLOSED-SHELL SYSTEMS 



An atomic many-body theory is size consistent if the properties calculated using it scales linearly as the number of 
the electrons. The diagonalization and the Bloch equation based methods are size consistent within a complete CSF 
space. But it is size inconsistent if the configuration space considered is incomplete. A consequence of incomplete 
cancellation of the unlinked terms, which scales nonlinearly to the number of electrons. The cancellation is complete 
when the CSF space is complete. To study atomic EDM which has important implications and small in magnitude 
it is preferable to use atomic-many body theory which is size-consistent even with an incomplete CSF space. It is 
also difficult to satisfy the condition of completeness for heavy atoms, which are important candidates for the atomic 
EDM experiments as the number of possible CSFs runs into millions for a moderate size orbital space. 



A. Size Consistency with Linked Diagram Theorem 



The wave-operator calculated using the Bloch equation within an incomplete CSF space is size-consistent if only 
the linked terms are retained |lq| . The incomplete cancellation of the unlinked terms is then avoided 



p = q[ v cs n cs p - Xcs pv cs n cs p 



linked 



Q v es n es p- Xes w 



(3.1) 
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where W = PV cs £l cs P. Redefine the wave-operator in terms of orders of excitation and consider only the single and 
double excitations. The wave-operator and correlation-operator are 



I + fles(l) + n es (2) = ^2 Qcs(m) and Xcs = ^ ^es(w), 



(3.2) 
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where m is the order of excitation. The closed-shell single and double excitation wave-operator diagrams are shown 
in Fig^. Similarly the V cs diagrams are shown in Fig^|. Unlike the single particle approach, the diagrams are used 
as representation of the physical effects and cannot be evaluated directly using the usual Goldstone rules as E-N 
partitioning is used. The W is the energy and hence a number for closed-shell sytems. From the definitions 



£l es (m),H 



Q(V es + Ksftes(l) + V cs n cs (2) - tt cs (m)W 



(3.3) 
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Let |$o) be the reference configuration and {I'E'q)} be the configuration space spanned by singly and doubly excited 
configurations. The wave-operators can then be expressed as 



IUl)=£k)(*oa£ and O e8 (2) = £ \*%){$. 
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(3.4) 



Where x r a and are the excitation amplitudes. For closed-shell systems with a single reference, W is just a number 
or closed diagrams. The term VL es W is therefore unlinked and does not contribute to the linked Bloch equation. The 
one-body wave-operator equation is 
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Similarly, the two-body wave-operator the equation is 
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Introduce H PTV as perturbation and redefine the perturbation Hamiltonian as V = VJ, S + H PTV . The corresponding 
wave-operators are 



0(1) = fies(l) + Oes,ed ro (l) and fi(2) = n es (2) + n csMm (2). 



(3.7) 



Wave-operators 6 s(l) and f2 C s(2) are same as before but fi e s,edm(l) and f2 eSje dm(2) connect to where 

is the configuration space opposite in parity to |$o)- Within the total configuration space r2 e s,edm(l) and 
^cs,cdm(2) can be represented as 
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The equations of O es ,edm(l) is 



(3.8) 
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The wave-operator defined by equation (3.9) has only one order of iJ PTV as terms of the form ($.;|iJ PTV \$j)x are 
excluded. The equation of fl PTV (2) is 
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The term ($ ah |_ff PTV |$o) does not contribute as the one-body interaction Hamiltonian i/ PTV cannot create double 
excitations. The equations (3.5)-( 3T0|) are the required wave-operator equations. 



B. Size Consistency with Connected Diagrams 



Consider the term ( $ 



€ ab ^a'/ x a' m (3-6) the diagrams of which are given in FigJ^. Among the diagrams (a) is 
linked but disconnected and remaining are connected. The disconnected diagrams in the wave-operator can introduce 
unlinked terms in the next iteration and hence a selection of only linked terms is required, this is difficult in terms 
of CSFs as all the contributions are combined. The other method of separating the wave-operator is in terms of 
connected cluster operators. It is easier to select connected terms, since it can be done without separating cluster 



operator into subcomponents. Consider the term (^hl^csl^a 



again, the disconnected contributions has a' and 



in and excitation is through the one-body part in V cs , avoiding these terms make the contributions from 



\^a') x a' connected. The third term in ( p,10[ ) still has disconnected terms but using the the same method 
these can be removed. After these modifications all the terms of the cluster equations are connected. To distinguish 
from the linked diagram excitation operators define the cluster-operator as T n , then 



A' 



T n = n(n) and T = ^ T n 



71=1 



Taking only the linear terms of one and two-particle cluster-operators, the wave-operator is 

n cs = 1 + T CS (1)+T CS (2) and W = V CS T. 
The wave-operator f2 0S is approximated by the linear cluster terms for the following reasons: 



(3.11) 



(3.12) 
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1. The correlation introduced by Tj 2 is very small compared to the contribution from T2, which represents a large 
part of the electron-electron correlation effect. 

2. Among the four-body cluster operators T% is the major contributor but in the present formalism this term can 
not be included as the CSF coupling is not in particle-hole form. 

3. Though T\ does not contribute significantly to the electron-electron correlation it is important since if PTV and 
the dipole operators are single-electron operators. 

The equation of the cluster-operator treating V cs as the perturbation is 



p = Qv cs n cs p - xosWP 



(3.13) 



In closed-shell systems XesW^ is always disconnected and do not contribute to the cluster equation 



T cs , Hq 



P = QV CB n es P 



(3.14) 



These are the CEPA-0 equations and do not include the EPV diagrams. The linked EPV diagrams should be avoided 
but unlinked EPV terms should be retained. After suitable transformations the unlinked EPV terms can be converted 
into connected terms |1S] , thus the cluster-operator equation is 



T cs , Hq 



P = QV cs n cs P 



EPV 



linked 
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Where the first term is EPV and second term is non-EPV. By rearranging 



EPV / % EPV 

QV.J-l-F) =-[Qx C sWP) 

linked V 



Then, the cluster-operator equations are 

T cs , Hq 

The one-particle cluster operator equation is 



EPO / ^ EPV 

P=[QV cs n cs P) -[QxesWP 



(3.15) 
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(3.18) 



The cluster amplitudes are denoted by T to distinguish from the one-particle cluster amplitudes represented by t. 
Here the calculation is using CSFs and T£ is the amplitude of the cluster operator which excites the reference CSF 
|$o) to the CSF 1$^). Similarly, the two-particle cluster amplitudes can be defined. 
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Similarly, the PT-violating cluster-operators T PTV are evaluated using the equations 
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Using the cluster-operators 

^es,edm — ^cs.cdm(l) ~t~ ^cs,cdm(2). 

The atomic EDM d a can be calculated using the operators as 
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(3.21) 



(3.22) 



(3.23) 



where Z) e ff = f2,t s -D^cs,cdm + ^cs cdm-^^ cs ^ s same as before except that the wave-operator is now in terms of connected 
clusters. The cluster equations are similar to the CEPA-2 equations (lj]]2l[], the explicit CEPA-0 equations are same 
as these equations without the EPV terms. 



IV. ANALYSIS OF THE CLUSTER EQUATIONS 



The cluster based formalism is not an order by order formalism but an iterative scheme where the Bloch equation 
is defined in orders of excitations rather than the orders of perturbation. It is possible to separate the contributions 
from various many-body effects cleanly using Moller-Plesset partitioning and an added advantage is the one-to-one 
correspondence with the diagrams using Goldstone evaluation rules. Such an approach has been used [^FJ to calculate 
parity non-conservation in atomic cesium to very high accuracy. The E-N partitioning mixes the contributions, but 
it has the advantage of capturing the static correlation effects very effectively. The CEPA-0 equations are exactly the 
linearized singles and doubles coupled cluster equations, the CEPA-2 equations includes a class of non-linear terms. 



A. The Singly Excited Amplitude Cluster Equation 

The diagrammatic representation of the principal terms in ( j3.1§| ) are shown in Fig.||. The contribution from each 
of the diagrams are as described: 

1. Diagram (a) contributes to the first term ($^|iJ cs |$o) of the cluster equation and is independent of any cluster 
amplitudes. It is an important term as the iteration proceeds from this term. 

2. The diagrams (b),(c),(d),(e) and (f) contribute to the second term (<& r a \H cs \<& r a ,)T£, and the V cs matrix element 
is coupled with the single excitation cluster amplitude. These start contributing from the second iteration, 
where the cluster amplitude is just the matrix element (<&^, | V^ s | 'I'o) in the first iteration. 

3. Diagrams (g),(h),(i) and (j) contribute to the third term {^ r a \V CB \^ r a , s b ,)T^,§ and couple the double excitation 
cluster amplitude with the V^ s matrix element. Similar to the second term, these diagrams start contributing 
from the second iteration. 

Consider the diagrams (c) and (e), though they resemble the Hartree-Fock potential scattering diagram these are 
very different. Consider the the bubble part of the diagram (c), it is summed over the occupied orbitals common to 
both the initial and the final CSFs in the matrix element of V cs . An example in Yb is if the initial and final CSFs 
are |3>£/f/) = 1 7s 2 ) and |3>£f) = |7s8s) respectively, the bubble part in (c) has all the occupied orbitals except the 6s 
orbital. This is because both the CSFs do not have 6s, where as in the Hartree-Fock scattering diagram the bubble 
should have contribution from all the occupied orbitals. A similar description is true of diagrams (d) and (f) too. 

The term (T a TT) EPV picks up the effect of the non-linear terms T es (l) 2 and T es (l)T es (2), which have T es (l) 
amplitudes. This implies that the wave-operator assumes the form 



n es = l+T 0S (l)+T es (2) 



T es (l)T es (l) + T os (l)T es (2) + T os (2)T es (2) 



EPV 



(4.1) 
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Terms which are not included in the single excitation cluster amplitude equation are 

T 1 EPO 

T CS (1)T CS (1) + T CS (1)T CS (2) + T CS (2)T CS (2) 

A later section describes the method to choose EPV terms from the renormalization part. 



(4.2) 



B. The Doubly Excited Cluster Amplitude Equation 



The diagrams of the principal terms of the double excitation cluster amplitude (3.19) are shown in FigJ|. The first 
term is similar to that of the single excitation cluster amplitude equation. The second term can be separated as 
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The first term on the right hand side has connected as well as disconnected terms, from which only the connected terms 
should be retained. The remaining terms are connected and hence linked too since the conditions a' ^ a, b; r' =/= ?*, s 
exclude the disconnected terms. Then 
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Similarly, the third term can be expanded to 



E (^|^s|^)^'=E E ( $ s|^e S |<^)^'+E 
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^ a'b' / * a'b' 



T ao / x ab 



where all the terms are connected. The triply and quadruply excited terms are excluded. Each of the diagrams has 
an exchange diagram too. Terms corresponding to each of the diagrams are: 

1. Diagram (a) correspond to the first term in the cluster equation and has no dependence on any of the cluster 
amplitude. 

2. (^ailVcsl^)^ contribute to diagrams (b) and (c). The final CSF in this term has a hole-particle pair in 
common with the cluster amplitude. 

3- {5 a ' a + <5a'h)( < I ) ab|V r C s| < I ) Q/)'7J contribute to diagram (c). Though the topology of the diagram is same as that of 
(^abl^esl^a)^! ^ is an EPO diagram. The V cs interaction changes the state of the particle and picks up a part 
of core-virtual correlation effect, which can be identified as core-polarization. 

4. (8 r >r + 5 r > s )(§ r a s b \V cs \$a,)T£' contribute to diagram (b). This also has similar topology with (^iK-sl^a)^ but 
is again an EPO diagram, where there is a change of the hole state and correspond to core-core correlation 
effect. 

5- (^allVcsl^ab^Tab' contribute to diagram (h) and (i). These are EPO diagrams where a hole-particle change to 
another hole-particle pair. These contribute to the core-virtual correlation effects. 

6. The term ( ( l ) ah|V r cs |$^ b ' s )72 b s contribute to diagrams (f) and (g). These are EPO diagrams where one of the 
particle states in T^ b s is excited to another particle state. This can also contribute to EPV diagrams of the first 
kind, if it is a hole-line EPV diagram then it will correspond to (h) and (i) and if it is particle line EPV then 
diagram (j). 
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7- ( < f > abl^ / cs| < i ) a6 S )^ai> S contribute to diagram (j). This is a double excitation where the particle states from the 
cluster amplitude T£ b s are excited to different particle states but the hole states remain intact. These terms 
capture the virtual- virtual correlation effects. 

8- (®ab\Vcs\&a'?t,)'l~a> S b contribute to diagram (d) and (e). These diagrams correspond to a change of the hole state 
and are EPO diagrams which capture the single-body hole-hole interaction component. This term can also 
contribute EPV diagrams, the hole line EPV diagram arising from this term is (k) and the particle line EPV 
diagrams are (h) and (i). 

9- (^abl^csl^l/b'^a'l' contribute to (j) and is a hole-hole correlation term. 

Thus terms in the cluster equation contribute to different physical effects. So far only the first three terms in the cluster 
equation have been considered. The last term in the doubly excited cluster amplitude (T^ b s W) Epy is a renormalization 
term. It picks up a set of terms non-linear in cluster amplitude 

T CS (I)T CS (2)+T CS (2)T CS (2). (4.4) 

These pick up a class of EPV terms which are non-linear in cluster amplitudes. Consider the expression of W 

W = PV cs (r cs (l)+T cs (2?jP. (4.5) 

The term T 0S (I)T CS (2) is picked up through PV CS T CS (1)P in W, which implies that (T a r fc ;s PV r cs T cs (l)) EPV can have 
one hole (particle) EPV line or a pair of hole-particle EPV lines. Whereas in the single excitation cluster amplitude 
equation, the contribution from T CS (1)T CS (2) is captured through the term (T a r PV CS T CS (2)) EPV in (T a r W) EPV . But the 
number of EPV hole-lines or EPV particle-lines are the same in both. In general the number of EPV hole-lines and 
particle-lines in (T es (n)PV es T es (m)P) EPY is limited by the V es if m,n > 2 and by the cluster amplitudes if m < 2 or 
n < 2. Although T es (l)T es (2) is included in the single as well as the double excitation cluster amplitudes the topology 
of the diagrammatic representations are different. Diagrams from (T^W) EPY has only a pair of hole-particle lines 
where as (T^ b W) EPY has two pairs of hole-particle lines. 



C. Selection of EPV Terms and Connected Terms 



The terms linear in cluster amplitude cannot violate Pauli exclusion principal. It is possible only when there are 
CSFs which are EPV, which is not possible. The EPV diagrams arise from renormalization terms, which are non- 
linear in cluster amplitude . The diagrams representing the renormalization terms are the cluster diagrams multiplied 
by the energy diagrams and hence unlinked but these can suitably be rearranged to yield connected diagrams. 

To select the EPV terms all the orbitals are tagged with labels which are prime numbers. The CSFs are also 
assigned a number which is the product of the prime numbers corresponding to the labels of the holes and particles 
of the CSFs. Consider a doubly excited CSF the prime numbers n a ,n b , n r and n s be the labels of the hole and 

particle states and N^ b their product, these five numbers identify the CSF. However, only three if the CSF is singly 
excited. To maintain consistency the remaining two indices are filled with another prime number not used in labeling 
the orbitals, let this number be Np but set the corresponding multiplying factor in n as unity. According to this 
scheme, the ground/reference state of Yb is identified by Np, N p , N p , A^andl. Similarly, labels are also given to the 
cluster amplitudes. 

The terms to retain from (T£W) EPY are those having CSFs in W identified by M which can be divided by one 
or more of the numbers identifying hole(particle) of the cluster amplitude T. The number of possible division is the 
number of common hole/particle lines between T and W. During the selection process division by N p should be 
discarded as this does not represent any hole of particle states and this is automatically achieved as the corresponding 
multiplying factor in n is unity. The advantage of this scheme is that it reduces the number of operations required in 
the selection process. 

The term (^allVcs]^,)^, of the doubly excited cluster amplitude equation has disconnected components if both 
the hole and particle states in the initial CSF are present in the final CSF. These are discarded and only the connected 
components are chosen. This can be implemented while calculating the matrix elements of V es . During the matrix 
element calculation the total number of hole states of the initial and final CSFs is calculated. If this is equal to three 
then these contribute to {^ r a b \V cs \^ r a i)T^, . In the next step the connected component is chosen by selecting the V es 
matrix element which has r' . Then the doubly excited cluster amplitude equation has only connected components. 
However the discarded components arc disconnected but linked. 
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V. THE CONFIGURATION SPACE CONSIDERED 



The configuration space is spanned by the CSFs constructed from the T^ N_1 orbitals. The CSFs are generated by 
single or double excitations from the occupied orbitals to the bound and the continuum virtual orbitals in all possible 
ways such that it yields the required final angular momentum. For the single reference MBPT, the reference CSF 
of Yb is |6s 2 ), which be referred as |$o)- Hence the occupied orbitals are (l-6)s,(2-5)p*, (2-5)p,(3-4)ci*, (3-4)d,4/* 
and 4/ respectively. The CSF space generated is not a complete active space but complete for the single and double 
excitations from the most important outer occupied orbitals within the orbital orbital space. The occupied-orbital 
shells of the configurations that has been considered are 



single excitation: 
double excitation: 



4/* 6 4/ 8 6s), 



4/* 5 6s ; 



4/* 5 6s 



4/* 6 4/ 8 
5p^4f 5 6s 
5s 1 5p* 1 6s 2 ) and 5s 1 5p 3 6s 2 



5p 3 4/* 5 6s 2 ), 5p* 3 4/ 5 6s 



4/ 7 6s 2 ), 5p* 1 6s 



5p 3 6s 2 y and 5s 1 6s 2 
4/ 7 6s), 4/* 4 6s 2 ), 4/ 6 6s 2 ^}, 4/* 5 4/ 7 6s 2 ) bp* 1 6s), 
5p 3 4/ 5 6s 2 \, 5s 1 6s 



5p 3 6s ), 



5s 1 4f* 5 6s 2 ), 5sHf 7 6s 



The remaining electrons are distrubuted among the virtual orbitals in all possible ways. From all the CSFs only the 
J = even parity and J = l odd parity CSFs are chosen. The number of non-relativistic CSFs generated are given in 
Table)! Although not included in the table, CSFs with excitations from 5s are also included in the CSF space. 

The total number of odd and even parity CSFs with bound virtual orbitals are 9930 and 17087 respectively. The 
modulus of the EN-partitioned energies of the CSFs-the diagonal Hamiltonian matrix elements-are as shown in the 
histograms Fig. 4. The two histograms are plotted such the the lowest \E\ is shifted to zero and the the range between 
the lowest and the highest are divided into ten units. The zero on \E\ axis are 14064.9531 and 14065.0068 hartrees 
for the even and odd parity CSFs respectively. Similarly, the highest \E\ are 14067.6720 and 14067.5996 hartrees 
respectively. From the histogram, the number of CSFs with low and high energies are less whereas the number of the 
configurations that can give the intermediate energy are large. As a result the perturbation series converges fast as 
only a few of the configurations are quite close to |$o) and the energy separation with the rest of the CSFs is quite 
large. 

Number of odd parity CSFs in the intermediate energy is more than the even parity as the odd parity configuration 
space can have many possible intermediate couplings to give J = l, to limit the number of CSFs within the memory 
limitations a selection of CSFs is done. The double excitations to d and / symmetries above the energy of the 
converged orbitals are not included. Another constraint on the choice of configurations is: there shouldn't be more 
than four open shells in the non-relativistic notation and eight in the relativistic form, choosing only singly and doubly 
excited configurations satisfies this condition for a closed-shell atom like Yb. This constraint is due to the angular 
co-efficient computation program. 



VI. RESULTS 



The atomic Yb EDM is calculated using the Bloch equation based MBPT and the size-consistent CEPA equations. 
As a part of the calculation the ground state energy is also computed. The results with different methods are given 
in the following sections. 



A. Bloch Equation Based MBPT 

1. Calculation ofQ es and Eo 
Using the wave-operator f2 os the ground state wave-function |^o) and energy Eq are 



^CS 




($0 


^CS^CS 





and E = ($o 



(6.1) 



The denominator in the expression of |^o) is the normalization factor. The first order energy correction of the ground 
state is zero Since the calculations are using the EN-partitioned Hamiltonian. 
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The even parity CSF space of the calculation is spanned by 9930 CSFs, of which core part of first 4435, 4436- 
9094 and remaining CSFs are |4/ 14 6s), |4/ 14 ), |4/ 13 6s 2 ) |4/ 13 6s) and |4/ 12 6s 2 ), and |5p 5 4/ 14 6s 2 ), |5p 5 4/ 14 6s) and 
|5p 5 4/ 13 ), and |5s5p 6 4/ 14 6s 2 ) respectively. Using this set of CSFs the ground state wave-function is 



0.9251 23 
-0.0600 60 
+0.0443 18 



6s' 



0.1172 17 



6p *' 



0.1169 21 



6s7s) + 0.0996 76 



5d" 



6p8p 



0.0568 61 



0.0442 15 



6p7pj - 0.0497 43 5cZ* 
6s8s\ - 0.0320 73 



0.0480 54 



6p* 7p* 



5d6d) 



+ 



Where only the ten important CSFs are given explicitly and the normalization constant is 1.0809 37. As expected, 
the most important CSFs in \^q) are doubly excited except for |6s7s), which is a singly excited CSF and does not 
interact very strongly with the ground state CSF |6s 2 ) but contributes significantly by correlation through other CSFs 
when the residual interaction is taken to higher orders. Values of Eq with increasing size of the even CSF space is 
given m Tablej|. CSFs are added to the calculation in sequence of excitations from the deeper core orbitals. As the 
size of the CSF space is increased more complicated many-body effects are included in the computation. 

An important quantity that can be extracted from Table|| is the change in Eq . Define the correlation energy AEq 
as the energy difference between the CSF energy of 1 6s 2 ) and the energy calculated using the CSFs in the even-parity 
CSF space. From the plot of AEq in Figfjja it is evident that the change in the ground state ASF energy is not 
uniform but in steps interrupted by regions of very minimal changes. Most significant changes of AEq occur while 
increasing the CSFs from 100 to 500, from 1000 to 2000 and from 6435 to 7435 CSFs respectively. These changes 
are not the combined effect of all the CSFs added but due to a few important ones. The largest change of AEq is 
while increasing from 1000 to 2000 CSFs which corresponds to contribution from the core configuration |4/ 14 ), that is 
double excitation from the 6s orbital shell. The remaining two are due to the core configurations |4/ 13 6s) and |5p 5 6s). 
Each contribute -0.0041 82, -0.0179 26 and -0.0045 07 hartrees respectively, the combined effect adds to 81.20% 
of the total correlation energy —0.0304 50 hartrees. From these it can be concluded that, the most important CSFs 
contributing to the correlation energy have core configurations |4/ 13 6s), |4/ 14 ) and |5p 5 6s). As to be expected the 
doubly excited CSFs are most important to capture the correlation effects and the low lying double excitations from 
6s orbital shell has the most significant contribution to correlation energy, it contributes 54.74% of the correlation 
energy. 

The plot (b) in Fig:^ indicates the need to include V es to high orders to capture the correlation effects accurately. 
From the graph the correlation effect due to two order of V es is —0.0436 033 hartrees and decreases in magnitude 
monotonically till fourth order to —0.0306 67hartrees but increases in the fourth order to —0.0343 38hartrees. This 
trend of oscillation about the final value of AEq continues till convergence. The cycle of the oscillation has a period of 
four orders, that is in four orders it goes to the same side of the final value of AEq and the amplitude of the oscillation 
decreases with each cycle. Over all, the value of AEq behaves like a damped oscillator with a cycle of four orders. 
If Eq is calculated by truncating the perturbation to the first few orders where the amplitude of oscillation is quite 
significant the value of A_Eo can be erroneous. 

The Fig.|| shows the trend of wave-operator convergence. The first graph Fig.^Ja is the value of the convergence 
criteria plotted against the order of perturbation and second graph Fig.||b is log 10 of the convergence criteria plotted 
against the order of perturbation. From the first graph it is evident that the convergence criteria begins with a small 
value but as shown in the second graph in terms of order of magnitude, the convergence is not so fast. The convergence 
is monotonic with very regular fluctuations. The wave-operator Q cs is stored in an order by order sequence. These 
are accessed as and when required during the calculation of r2 es ,edm- 



2. Calculation of Sl OSjCt i m and d a 

The wave-operator f2 C s,edm is calculated introducing the interaction Hamiltonian 7J PTV and adding opposite parity 
CSFs to the CSF space. Once H PTV is applied to the wave-operator f2 cs , it maps onto the odd-parity component of 
the CSF space and can never be mapped back to the even-parity space as H PTV is treated to first order only. This 
is followed by a sequence of residual coulomb interaction V^ s , which accounts for the correlation effects within the 
odd-parity sub-space. In sum total it is a sequence of perturbations applied to the ground state CSF |6s 2 ), where 
H rTV is sandwiched between all possible arrangements of V cs . After the wave-operator r2 cs .edm is calculated, the mixed 
parity ground state wave- function \^q) can be written in terms of the wave-operators Q cs and f2 C s.cdm as 

*o) - 



#0 



a 



es. cam 



(6.2) 
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The value of d a can be calculated using the expression 



d„ 



D 



2(*r 



(6.3) 



Choosing the odd parity CSF as mentioned in Sec. V the correction to the ground state |\I>o) from the opposite parity 
sub-space due to H PTV is 



#corr)=-A( - 55.1403 73 
-9.5064 23 
+5.5464 22 



6s6p*J - 17.3681 84 
hp* 7s^ - 9.2841 00 
6s9p*\ + 5.5189 58 



6s7p*^> + 10.8231 86 
5s6p*) + 7.8322 72 
6p5d*) + 5.1246 77 



6s6p 

6s8p* ) + 7.1950 5(1 

5p* 8s 



6p* 5d ■■ 



Where A — \Z2Ct&nG f and only the first ten important CSFs are listed. The coefficients are much larger than 
unity but these should be scaled by the parameter V2G F . The product of the coupling constant Ct and nuclear spin 
ctjv is retained as a parameter and Ct can be estimated after combining with the experimental results. The above 
expression for l^^n.) shows that: 

1. Like in the lowest order single-particle calculation, the coefficients of the CSFs \6snpt) flip sign for n>8. 

2. Most of the important CSFs are singly excited with respect to the ground state CSF |6s 2 ), which is to be 
expected since H FTV is a single particle interaction Hamiltonian. 

3. Singly excited configurations like |6s6p) can contribute through three possible many-body routes. First H PTV 
excites |6s 2 ) to \6snp*) then a sequence of V es connects it to |6s6p), second a sequence of V es connects 1 6s 2 ) to 
\6p* 6p) and iZp TV takes it to \6s6p) and third a sequence of V cs takes |6s 2 ) to a CSF which connects to 
|<&j) via 7J PTV and another sequence of V cs connects it to |6s6p). Although |6s6p) cannot connect directly to the 
ground state |6s 2 ) as 6p is close to zero within the nucleus, it is the third most important CSF which contributes 
to l^corr)- This demonstrates the importance of many-body effects. 



4. The two most important doubly excited odd-parity CSFs for the evaluation of EDM are 1 6p+=5t^ and \6p5d*) . More 
interesting is the second as both the virtual orbitals involved cannot contribute to the H rTV matrix elements. 
Among the possible many-body routes which can contribute to the co-efficient of \6p5d*} one possibility is 
through the deeper occupied orbitals 5p* and 5s, which would contribute to core polarization effects. 

The value of d a is calculated using different sets of CSFs, where the number of either the even or odd CSFs are fixed 
to the maximum allowed and then increase the number of CSFs in the opposite parity CSF space. The results of such 
a sequence of calculations is given in the Table III . If the previous sequence of calculation shows the importance of 



occupied orbitals in the whole CSF space, these two sequences demonstrates the significance of the occupied orbitals 
in CSF sub-space of each parity. 

Like the correlation energy AEq there is a significant change of d a when CSFs with double excitations from 6s are 
included. To appreciate the change better the values of d a in the two sequence are plotted in Fig.||. Consider the 
sequence where the number of even-parity CSFs is fixed, d a increase with the number of the odd-parity CSFs, which 
implies that there are no appreciable cancellations due to the odd parity CSFs added. In the second sequence where 
the number of odd-parity CSFs is fixed, d a decreases as the number of even-parity CSFs is increased. 

In both the sequence there is a significant change in d a when CSFs with double excitation from 6s orbital are added. 
But, the changes are different in sign, in the even CSF space the inclusion of CSFs with double excitation from 6s 
increases the value of d a where as in the odd-parity CSFs it decreases. Consider the expression for EDM it can be 
expanded as 



2 $ 



'-cs.cdm 



$0 ) + H( $0 



cs.cdm 



$0 



Which can be rewritten in terms of CSF coefficients as 



,(es) 



Da 



cs.cdm 



2 E H + E(*» 



(es) 



DC 



(cs,cdm) 



(6.4) 



(6.5) 
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Where the definitions of all the quantities are the same as defined in paper-I. Within the whole CSF space, the 
contribution from the first term is 5.4394 39 and the contribution from the second term is —0.6417 31, which is just 
11.80% of the first term. Five most important configurations in the second term from the even-parity sub-space are 
|6p* 2 ), |6p 2 ), |6p*7p*), \5d* 2 ) and \6p7p) and their contributions are -0.8119 50, 0.2065 15, -0.1001 58, 0.0738 55 and 
0.0599 84 respectively, where d a is in units of Ct<Jn x 10~ 12 eao- All these are doubly excited CSFs and mixes with 
the ground state CSF significantly but the singly excited CSF |6s7s) which is the third important CSF of |^o) does 
not contribute strongly. In addition there is shift in the sequence of the important CSFs compared to the sequence of 
contribution to |\&o)j this is due to difference of dipole and H PTV coupling strengths between different CSFs. 

Within the whole CSF space considered the value of d a is 4.4438 58. In absolute terms this is 0.4804 56 less 
compared to the lowest order result of 4.9243 136 calculated in an earlier section. A major contribution to this 
difference is the many-body effects, this is because the direct contributions from the CSFs added to the configuration 
space is small. Which implies that the contribution from the many-body effects is just 10.81% of the total value 
and the change is negative. An added advantage of the order by order approach is that the contribution to d a can 
be calculated in terms of the order of residual coulomb interaction. Earlier while calculating the ground state ASF 
energy it was shown how a truncation in the order of V cs perturbation can give an inaccurate value of Eq. 



B. Cluster Based Formulations 



The CEPA-0 wave-operator equations of Q es and £l e s,edm are identical to the linearized coupled-cluster equations. 



The atomic Yb EDM d a calculated using CEPA-0 using the same sequence of CSFs as in Table III are given in Table 



IV and the following can be inferred: 



1. The results from CEPA-0 do not differ significantly from the MBPT results for the singly excited CSFs with core 
configurations 1 6s) and |4/ 13 6s 2 ). This is because the non-linear terms do not contribute since the configuration 
space is spanned by singly excited CSFs. 

2. A significant difference from the MBPT results is expected when doubly excited CSFs are included in the CSF 
space. This is observed when the CSFs having core configuration |4/ 13 6s 1 ) are included. The change of d a 
calculated using MBPT by fixing the number of odd parity CSFs to 17087 is -0.5599 84, whereas CEPA-0 
reduces the change to —0.1091 50. 

3. Within the whole CSF space considered, the value of d a calculated using CEPA-0 is 5.9421 36, which is 25.21% 
larger than the MBPT result. The difference gives a rough estimate of the contribution from the non-linear 
terms in the cluster equation. It is a rough estimate as the MBPT calculation also includes the contribution 
from the size-inconsistent terms, which are unphysical and are not included in the CEPA formalism. 

The CEPA-0 calculation does not include size-consistent non-linear terms of f2 es and £l e s,edm- The results using the 
CEPA-2 equations are tabulated in Table [v| To make comparisons convenient the calculation is done with the same 
choice of CSF sequence as before. Comparing with the results from the calculations using MBPT and CEPA-0 give 
the following: 

1. The calculation within the CSF space of singly excited even parity CSFs with core configurations 1 6s) and 
|4/ 13 6s) and all the odd parity CSFs is identical with the result calculated using CEPA-0. This implies that 
the contribution from the non-linear EPV terms with T es (l) from the even parity subspace is negligible. But, 
it is different from the MBPT result of the same CSF space. And the difference of 0.2067 67 is due to the EPO 
non-linear terms, which are included in the MBPT. 

2. Consider the calculation within the CSF space consisting of all the even parity CSFs and the singly excited odd 
parity CSFs having core configurations |6s) and |4/ 13 6s). The CEPA-2 and MBPT results differ by less than 
2%. But, the result from CEPA-0 differs from both by more than 12%. This shows that the contribution from 
non-linear terms with Tptv(1) is not negligible or the contribution from the size-consistent is significant, with 
the present results it is difficult to distinguish between the two. 

3. The d a calculated within the whole CSF space considered is 4.5065 25 and is larger than the MBPT result by 
1.4% and less than the CEPA-0 result by 31.86%. This implies that the contribution of non-linear terms to d a 
is very important. 

The effect of the EPO unlinked terms in f2 e s,edm can be estimated by calculating T es ,edm with the EPO rcnormalization 



terms included. The results are given in Table VI. The final value with the full CSF space is 4.2446 91, which 
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is suppressed by 5.8% compared to the CEPA-2 result. This difference is due to the EPO contribution to the 
renormalization term in fi e s,cdm- Another property which can be compared to gain an insight on the contribution of 
the size-inconsistent terms is the energy of the ground state E . The value of E calculated using the wave-operator 
fi es derived here has no contribution from the size inconsistent terms but it excludes some of the less important 
size-consistent terms. The difference in the value Eq calculated using the Bloch-equation and Q cs derived from the 
cluster equation gives the contribution from the size inconsistent terms. Like in d a the approximation is that the 
contribution form the EPO size consistent terms non-linear in cluster amplitudes is very small, then the difference in 
the result can in principle be accounted to the size inconsistent terms. 



VII. CONCLUSION 



Comparing the results from different methods, it is clear that a size-consistent theory is preferable for high accuracy 
computation of atomic EDM and the contribution from the non- linear terms in the cluster amplitude is also important. 
With the present calculation the difference between the MBPT and CEPA-2 result cannot be accounted distinctly to 
non-linear terms or the size-inconsistent terms. A better comparison can be made after including all the non-linear 
terms in the cluster amplitudes. As the orbital space is made larger, the size of the CSF space grows very large. This 
puts a limitation during the calculation as the memory requirement increases, it is manageable if the calculation is 
done at the level of single particle Moller-Plcssct partitioning of the atomic Hamiltonian. At the single particle level 
the cluster equation with the residual coulomb interaction Hamiltonian reduces to the method used used by Blundell 
|p(if and his collaborators. The inclusion of terms non-linear in cluster amplitudes is also relatively easy as compared 
to the EN-partitioning. To check the quality of the wave functions, experimentally known quantities like hyperfme 
constants and excitation energies needs to be calculated. These will be reported in our later papers. 
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(a) (b) 
FIG. 1. Diagrams for the wave-operators (a)57 es (l) and (b) r2 es (2) . 




(n) (o) O < m > 

FIG. 2. The diagrams for the residual Coulomb interaction V ea - 




(a) (b) (c) 

FIG. 3. The diagrammatic representations of the term {&ab\Vcs\&ai}Xa, 
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The energy of the even parity configurations. 



The energy of the odd-parity configurations 



(a) (b) 
FIG. 6. Histogram of \E\ for the even parity CSFs. 
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(a) (b) 
FIG. 7. The change in energy due to many-body effects introduced by the configurations. 
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(a) (b) 
FIG. 9. The value of d a 



TABLE I. The number of the CSFs with different occupied configurations. 



Sl.no 


Occupied part 


Configurations 


Sl.no 


Occupied part 


Configurations 






Even 


Odd 






Even 


Odd 


1 

3 


|6s> 

|4/ 13 6s 2 > 
|4/ 12 6s 2 > 


6 
12 


12 
18 


2 
4 


|4/ i4 > 


147 
1224 


287 
3618 


5 


3045 


7739 


6 


|5p 5 4/ 14 6s 2 } 
|5p 5 4/ 13 6s 2 } 


12 


30 


7 




1044 


1668 


8 


3604 


2394 
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TABLE II. The energy of the ground state ASF with increasing CSF-space size. 



SI. no 


No oi Cbr 


Energy 


SI. no 


No oi Cor 


Energy 


i 
i 


i nn 

1UU 




o 
A 


ouu 




3 


1000 


-14067.6756 60 


3 


2000 


-14067.6942 64 


5 


3000 


-14067.6942 68 


6 


4000 


-14067.6949 91 


7 


4435 


-14067.6949 91 


8 


5435 


-14067.6974 13 


9 


6435 


-14067.6974 21 


10 


7435 


-14067.7019 28 


11 


8435 


-14067.7040 97 


12 


9094 


-14067.7040 97 


13 


9594 


-14067.7042 26 


14 


9930 


-14067.7042 26 



TABLE III. Values of d a for different number of even and odd parity configurations. 



SI. no 


No of CSFs 


d a 


SI. no 


No of CSFs 


d a 


Even 


Odd 


Even 


Odd 


1 


9930 


12 


3.3654 83 


10 


7 


17087 


8.2430 89 


2 


9930 


30 


3.3605 54 


11 


19 


17087 


8.2430 89 


3 


9930 


3648 


3.4095 46 


12 


1243 


17087 


7.6831 05 


4 


9930 


3846 


4.3895 21 


13 


1390 


17087 


4.2645 44 


5 


9930 


11676 


4.3895 43 


14 


4435 


17087 


4.2118 95 


6 


9930 


14070 


4.3916 11 


15 


4448 


17087 


4.1837 51 


7 


9930 


14100 


4.3760 30 


16 


6713 


17087 


4.5282 75 


8 


9930 


15768 


4.4424 65 


17 


9094 


17087 


4.4333 81 


9 


9930 


17087 


4.4438 58 
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TABLE IV. Value of d a computed using the CEPA-0 formalism. 



SI. no 


No of CSFs 


d a 


SI. no 


No of CSFs 


d a 


Even 


Odd 


Even 


Odd 


1 


9930 


12 


3.8250 23 


10 


7 


17087 


8.4499 64 


2 


9930 


30 


3.8426 96 


11 


19 


17087 


8.4499 64 


3 


9930 


3648 


3.9252 41 


12 


1243 


17087 


8.3408 14 


4 


9930 


3846 


5.7783 09 


13 


1390 


17087 


5.5126 26 


5 


9930 


11676 


5.7783 88 


14 


4435 


17087 


5.5091 16 


6 


9930 


14070 


5.7820 42 


15 


4448 


17087 


5.3769 11 


7 


9930 


14100 


5.8370 85 


16 


6713 


17087 


5.6132 71 


8 


9930 


15768 


5.9437 56 


17 


9094 


17087 


5.9275 26 


9 


9930 


17087 


5.9421 36 











TABLE V. Table V: Values of d a computed with CEPA-2. 



SI. no 


No of CSFs 


d a 


SI. no 


No of CSFs 


d a 


Even 


Odd 


Even 


Odd 


1 


9930 


12 


3.3130 75 


10 


7 


17087 


8.4499 64 


2 


9930 


30 


3.3133 86 


11 


19 


17087 


8.4499 64 


3 


9930 


3648 


3.3684 39 


12 


1243 


17087 


7.8528 23 


4 


9930 


3846 


4.4589 69 


13 


1390 


17087 


4.3791 81 


5 


9930 


11676 


4.4589 33 


14 


4435 


17087 


4.3758 71 


6 


9930 


14070 


4.4608 15 


15 


4448 


17087 


4.3558 36 


7 


9930 


14100 


4.4581 78 


16 


6713 


17087 


4.4547 38 


8 


9930 


15768 


4.5149 07 


17 


9094 


17087 


4.4981 25 


9 


9930 


17087 


4.5065 25 
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TABLE VI. Values of d a computed with the EPO components included in the cluster amplitudes T. 



SI. no 


No of CSFs 


d a 


bl. no 


No of CSFs 


d a 


Even 


Odd 


Even 


Odd 


1 


9930 


12 


3.1704 91 


10 


7 


17087 


8.4499 64 


2 


yyou 


OA 

oU 


3.1682 67 


11 


ty 


1 (Uo 1 


O A A C\C\ /"> A 

8.4499 64 


Q 
O 


9930 


3648 


'X 91 Oft ^ 


1 9 


1243 


17087 


7 8^91 8/1 


4 


9930 


3846 


4.2130 96 


13 


1390 


17087 


4.3783 45 


5 


9930 


11676 


4.2130 54 


14 


4435 


17087 


4.3042 71 


6 


9930 


14070 


4.2147 05 


15 


4448 


17087 


4.2332 39 


7 


9930 


14100 


4.2039 41 


16 


6713 


17087 


4.2802 30 


8 


9930 


15768 


4.2539 59 


17 


9094 


17087 


4.2375 23 


9 


9930 


17087 


4.2446 91 
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